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Представлена модель учета деформаций в узлах трения типа «колесо-опора». Модель позволяет оценить 
влияние динамики системы управления на степень деформации контактирующих поверхностей И наоборот. 
Предложено описание возникновения ВОЛНОВЫХ деформаций рельс в железнодорожном транспорте как 
следствие действия циркуляционных сил. 
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Введение. Оценка динамики сложных современных систем невозможна без учета процессов, 
возникающих в результате действия сил трения. При этом для многих систем трение является не 
каким-то побочным явлением, от которого по возможности избавляются, а процессом, которым 
можно и нужно управлять. Управление невозможно без разработки адекватной математической 
модели, отражающей основные закономерности и характеристики процесса, которым предстоит 
управлять. 

Процесс формирования и разрушения фрикционного контакта в узлах трения зависит от 
динамики движения самой системы, которая, как известно, определяется фазовыми траектория- 
ми. Оценка влияния фрикционного контакта на динамику систем управления носит эмпирический 
характер. Так, к примеру, в железнодорожном транспорте формирование силы тяги локомотива в 
узле «колесо — рельс» описывается нелинейной функцией, полученной усреднением статистиче- 
ских данных по динамике локомотива на заданном отрезке пути [1]. Такое представление о силе 
трения носит статический характер, не учитывающий динамику формирования и разрушения 
фрикционного контакта. Использование статической теории в моделях движения потребовала 
разделить тангенциальное усилие, возникающее при фрикционном контакте на силы трения: ка- 
чения, сцепления, скольжения. Такая дифференциация оправдала введение коэффициентов тре- 
ния (усредненных величин), позволивших, в свою очередь, существенно упростить расчеты. Но 
модели систем управления, построенные на этих приближенных оценках, не позволяют объяснить 
многие процессы, имеющие место в системах с трением, такие, к примеру, как фрикционные ав- 
токолебания или появление волновых деформаций рельсов. 

Деформационные изменения в зоне контакта можно учесть за счет введения соответст- 

вующей фазовой координаты, такое изменение повысит порядок системы дифференциальных 
уравнений, описывающих движение системы управления, но при этом позволит связать динамику 
системы и динамику контакта. 
Анализ формирования силы тяги в контакте колеса с поверхностью дороги. Формиро- 
вание вектора силы тяги в колесном транспорте, в том числе и железнодорожном, происходит в 
зоне контакта колеса с поверхностью дороги. В железнодорожном транспорте определяют зави- 
симость силы тяги от величины продольного проскальзывания колеса относительно поверхности 
рельса [1]. Продольное проскальзывание (крип) является разностью между скоростью движения 
поверхности колеса и скоростью перемещения центра масс колеса относительно поверхности до- 
роги. В условиях чистого качения (по инерции) эти скорости совпадают. Зависимость между тан- 
генциальным усилием и крипом представлена на рис.1 [1]. 
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Рис.1. Зависимость тангенциальных усилий от проскальзывания 


При анализе данной зависимости возникает закономерный вопрос, откуда берется про- 
скальзывание, чем оно обусловлено? Проскальзывание - явление, связанное с упругими и упруго- 
пластинчатыми деформациями в материале дороги (рельса) [2]. Так, анализируя зависимость, 
представленную на рис.1, можно увидеть, что часть характеристики до точки насыщения являет- 
ся зоной упругих деформаций, точка насыщения характеризует модуль упругости, а часть харак- 
теристики за точкой насыщения является зоной упруго-пластинчатой деформации. 

Возникновение упруго-пластинчатых деформаций в материале рельса под воздействием 
тангенциальной силы показано на рис.2. 
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Рис.2. Возникновение продольного проскальзывания деформаций: 
а - рельса в зоне контакта колеса и рельса; 6 -— пути под действием тангенциальной силы. 


На представленном рисунке рассмотрен крайний случай, когда тангенциальное усилие 
превышает модуль упругости. В большинстве случаев, когда движение происходит в зоне дейст- 
вия силы трения сцепления (ЁР‚„„. до точки насыщения, см. рис.1), тангенциальное усилие не пре- 
вышает модуль упругости. В материале рельса под действием тангенциальной силы происходит, с 
одной стороны, растяжение, а с другой стороны, сжатие материала. Отметим, что в этом процессе 
участвует только верхняя часть рельса, при этом глубина слоя участвующего в сдвиговых дефор- 
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мациях определяется степенью деформации колеса и рельса в зоне контакта под воздействием 
нормальной силы. 

Для построения модели тангенциального движения центра масс колеса, необходимо 
учесть влияние на параметры такого движения динамики деформационных перемещений центра 
масс под действием нормальной силы. 

Модель вертикальных перемещений без учета относительного движения контакти- 
рующих поверхностей. Рассмотрим систему, оказывающую влияние на контакт, силу и связан- 
ную с ними систему координат, на примере взаимодействия колеса с поверхностью дороги. 


ги ы 





Рис.3. Система координат и силы, связанная с центром масс 


На рис.3 представлено колесо (радиусом №) и рельс (опора). Силы и моменты, приведен- 
ные к центру масс колеса: Р — сила тяжести, № -— сила реакции опоры, М — внешний вращающий 
момент, Р.„и- сила, вызывающая линейные перемещения центра масс колеса (вдоль оси х), зави- 
сящая от момента М. 

Под воздействием силы тяжести Р центр масс колеса начнет движение вниз (положитель- 
ное направление оси у), деформируя колесо и опору, при этом сила, действующая навстречу это- 
му движению Л! (сила реакции опоры), начнет нарастать, по мере увеличения деформации колеса 
и опоры (координата у) до тех пор, пока не уравновесит действие силы тяжести. Нам не известны 
характеристики силы реакции опоры, но известно, что она изменяется от нуля до значения равно- 
го силе тяжести и это изменение носит характер, близкий к апериодическому [3], следовательно, 
будет справедливо следующее равенство: 

2 
р (1) 
4 Ч 

В выражении (1) т- приведенная к контакту масса колеса в кг, Й, с — характеристики вяз- 
кости и жесткости деформируемого слоя материала рельса и колеса в Н.с/м, Н/м. 

Если к колесу приложен внешний вращающий момент (см. рис.3), то возникает враща- 
тельное движение и связанное с ним, через трение в контакте колеса и опоры, линейное пере- 
мещение центра масс колеса. Не проводя анализ силы трения, оценим влияние этих двух видов 
движения на деформацию контактирующих поверхностей. Данные, полученные опытным путем 
[3, 4], показывают, что зависимость деформационных изменений от скорости вращения колеса и 
линейного перемещения имеет вид близкий к апериодическому, при этом угловая скорость колеса 
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влияет на деформацию зоны контакта колеса, 

а линейная скорость влияет на деформацию а 7 

зоны контакта опоры. ий 
Рассмотрим это явление на примере, Сцепление ; 


приведенном на рис.4 [1]. Здесь представле- Микроскольжение, 
> 


но изменение пятна контакта на рельсе при 
И. 





увеличении относительной скорости сколь- Е. .-:1.1-: , 

жения, движение системы происходит на вос- Микроскольжение / = 
ходящем участке кривой. С увеличением ско- Сцепление „” 1 
рости происходит как бы всплывание колеса / 

и пятно контакта существенно уменьшается, Микроскольжение/ 


при некотором критическом значении сцеп- 
ление теряется вообще и остается только 
скольжение. 

В случае значительной деформации 
опоры и колеса, при наличии линейного перемещения центра масс, колесо движется не только 
тангенциально (горизонтально), но нормально (вертикально). Вертикальное перемещение центра 
масс колеса связано с тем, что при движении колеса происходит замена контакта деформирован- 
ных поверхностей на недеформированные части. Используя предложенный ранее подход по нор- 
мировке действующих в зоне контакта сил относительно нормальной нагрузки, вместо выражения 
(1) получим следующую зависимость: 

ах 2 
Ра 


Скольжение/ 


Рис.4. Переход от сцепления к скольжению 





ау 
+ й, + 2 
ЧЕ? ПЕ (2) 
до % 


где с, «4 соответственно угловая и линейная скорости колеса в рад/с, м/с; Аз, Ал — коэффи- 


циенты, зависящие от свойств материалов, из которых состоят колесо и опора. 
Предлагаемый подход к моделированию влияния угловой и тангенциальной скорости 
движения на вертикальные деформационные движения центра масс колеса, похож на действие 
выталкивающей гидродинамической силы. Следуя таким рассуждениям, центр масс колеса при 
увеличении скорости движения (угловой и тангенциальной) начинает как бы всплывать, при этом 
уменьшается деформация контактирующих поверхностей и как следствие этого, уменьшается си- 





ла трения. 
Раскроем в правой части выражения (2) скобки и приведем подобные, получим: 
аа ах 2 
чи. >) а а 
Ре ‘течь че. (3) 
ИИА 1 
Ах 

Для упрощения выражения (3) учтем тот факт, что тя Иа С -) при этом из геометрии 


_ 
известно, что длина дуги окружности линейно зависит от угла — 1= 180’ а эффект относитель- 


ного проскальзывания контактирующих поверхностей колеса и дороги, в условиях стационарно- 
сти тяговой характеристики, можно линеаризовать коэффициентом Аз, значение которого нахо- 
дится между нулем и единицей. Исходя из этого, представим выражение (3) следующим образом: 


4”у ау 
т + 4 
4? НЕ (9) 





о 
Ре": =т 
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где К; = К, т © в 

180 41 
Связь вертикальных деформаций контактирующих поверхностей с тангенциальными 
силами в зоне контакта. Принято разделять действующую в зоне контакта тангенциальную 
силу, возникающую в условиях вращательного движения колеса, как совокупность следующих 
сил: 

— сила трения сцепления (покоя) Ру. сц.; 

— сила трения скольжения Р;рк; 

— сила трения качения Р;рх.. 

Сила трения сцепления характеризует участок кривой, представленной на рис.1, до точки 
насыщения, сила трения скольжения характеризуется падающим участком этой же кривой. Сила 
трения качения, имеет различную интерпретацию, но наиболее общим будет высказывание, что 
трение качения - это потери на продольное движение центра масс при условии, что тангенци- 
альная сила будет находиться на участке до точки насыщения (см. рис.1) [1]. 

Учет нормального деформационного движения центра масс колеса позволяет легко ин- 
терпретировать возникновение потерь при продольном движении центра масс колеса как потери 
тангенциальной силы, пропорциональные текущему значению нормальной деформации «у» 


(рис.5). 


у 


Рис.5. К объяснению потерь тангенциальной силы 


Для того чтобы центр масс колеса начал перемещение вдоль оси х необходимо перенести 
нормальную нагрузку из точки (6) в точку (а), что требует от тангенциальной силы провести ра- 
боту, обратную работе нормальной силы по деформации колеса и опоры в точке (6). Иными сло- 
вами тангенциальная сила должна поднять центр масс колеса вертикально вверх на расстояние у 
(текущее значение деформационных изменений колеса и опоры). 

В предлагаемом подходе к описанию силы трения качения (потерь на продольное движе- 
ние) нет необходимости разделения тангенциальной силы на силы трения покоя и скольжения. 
Так как в обоих случаях будут присутствовать «анти — деформационные» потери пропорциональ- 
ные у. Дифференциальное уравнение, описывающее продольное движение центра масс, будет 
иметь следующий вид: 

Я ах 
а" ар, 





(5) 


тяг 


Ах 
Е — УуР51еи(—) = т 
у Е 


где, К 


тяг 


р Ах 

- тангенциальная сила, приложенная к центру масс колеса Н, УР51{2п —. ) - потери 
ГА 

тангенциальной силы на преодоление деформационных изменений контактирующих по- 


. 4х 
верхностей Н, 512 п о - нелинейная функция, связывающая потери на трение качения с 
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продольным движением центра масс колеса, т ой - параметры, характеризующие инер- 


ционные и вязкие свойства продольного движения центра масс колеса. 
Тангенциальная сила (Р’..) зависит от текущего значения нормальной деформации кон- 


тяг 


тактирующих поверхностей колеса и дороги «у», так как, чем больше деформация, тем сильнее 
связь контактирующих поверхностей друг с другом. Текущее значение - «у» удобно нормировать 
относительно максимального, и полученное нормированное значение в дальнейшем использовать 
как весовой коэффициент при тангенциальной силе: 


о о (6) 
У тах 


где у„., определяется отношением нормальной силы к значению ст. 


Л . ах 
Применение нелинейной функции связи И) в уравнении (6) обусловлено тем, что 
1 
ах 
потери на трение качения присутствуют только в случае наличия движения, т.е. если ар #0. 
1 


: . ах 
Вместе с тем, функция 512" крайне неудобна, так как в точке т = 0 она имеет разрыв. Удоб- 
ГА 


но заменить функцию 51епи на гиперболический тангенс: 


Кх — Ах 
е” -е 
Ш (х) т НУ (7) 
ее” +е 
где А — коэффициент характеризующий степень близости функций 5еи (х) и Ш(х), 
тИй(х)] =9$1еи(х). (8) 
Гиперболический тангенс -— гладкая функция, непрерывная на интервале 


хе (—0 ь со) ‚ При этом степень пологости функции вблизи точки ноль определяется значением 


К. 
С учетом предложений и замечаний (6-8) выражение (5) примет следующий вид: 


2 

РЖИ Г о - би 
Е таг.тек. ЧЕ АЕ? "р 

Таким образом, выражение (9) описывает линейное продольное движение центра масс 
колеса в условиях, когда вертикальное движение центра масс под действием нормальной силы 
описывается выражениями (2), (3). 
Возникновение циркуляционных сил в зоне контакта колеса и опоры. На рис.2, 6 пред- 
ставлен участок железнодорожного пути с продольными деформациями, которые носят колеба- 
тельный характер. Рассматриваемая модель учета влияния деформационного сближения контак- 
тирующих поверхностей колеса и опоры на динамику движения центра масс колеса позволяет 
достаточно просто (для случая одной массовой системы) объяснить это явление как следствие 
действия циркуляционных сил. Рассмотрим систему уравнений, описывающих нормальное и тан- 
генциальное движение центра масс колеса в рамках предлагаемого подхода: 


(9) 








2. 
ее. ору И те 
а таг.тек. ДЕ ДЕ пр ДЕ (10) 
Р бани, 6) 4?у р у | 
е ь 2 + —+Сс 
ЧЕ? В" 


Реакция системы (10) на изменение координаты «у» вызовет возникновение циркуляци- 
онной составляющей вектора силы и как, следствие, процессионное движение центра масс коле- 
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са. Докажем это, проведя анализ матрицы жесткости системы (10). Для начала проведем не- 
сколько замен в (10), позволяющих указать связь сил, действующих на систему координат движе- 
ния центра масс, получим следующую систему уравнений: 


ах ах 
о 








пр Е 
4? 4 нА (0 
Е, =т р Ес 
2 а? 1 а 1У 
ах Сб 
де К (у) = Е... = УР) ,а Е, = Ре ыы 
Матрица жесткости системы, описываемой уравнениями (11), будет иметь следующий вид: 
0 ОР, 
С= ду |, (12) 
0 С 
ОЕ. 1 ах 
где “= ( И „тек. — Рй (—)) . 
СО 4 


Матрица (12) несимметрична, путем простых преобразований такая матрица превращает- 
ся в сумму симметричной и косо-симметричной матриц: 
ОЕ ОЕ 
0 0.5" 0 0.5—" 


д д 
се и Ни (13) 


0.5 ой С —0.5 — 0 
ду бу 

В выражении (13) первая матрица симметрична относительно главной диагонали и харак- 
теризует вертикальную и горизонтальную составляющую сил, действующих в контакте колеса и 
рельса, а вторая матрица косо-симметрична и характеризует циркуляционные силы, которые, в 
свою очередь, вызывают волнистые деформации в материале рельса (см. рис.2, 6) 
Заключение. О существовании деформационных изменений контактирующих поверхностей ко- 
леса и опоры известно давно, но учет этих изменений крайне сложная и неоднозначная задача, 
решаемая многими авторами ранее [3, 4]. Отличие предлагаемой модели от известных [3, 4] со- 
стоит в том, что она позволяет связать деформационные изменения в контактирующих поверхно- 
стях с основными параметрами движения системы не на каком-то этапе контактного взаимодейст- 
вия, к примеру, в режиме скольжения, а во всех режимах. Такой подход позволит моделировать 
движение систем, имеющих фрикционный контакт типа: «колесо-опора» во всем диапазоне изме- 
нений фазовых координат. 
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